TRANSPORT AND CONCENTRATION PROBLEMS 
WITH INTERACTION EFFECTS 



F. SANTAMBROGIO 

Abstract. After presenting an overview about variational prob- 
lems on probability measures for functionals involving transport 
costs and extra terms encouraging or discouraging concentration, 
we look for optimality conditions, regularity properties and ex- 
plicit computations in the case where Wasserstein distances and 
interaction energies are considered. 



1. Introduction 

One of the aim of the paper is to give a short overview of possible 
variational problems involving transport costs between distributions of 
mass and their concentrations. The general problem we are interested 
in is 

min v) := T(u, v) + F(fi) + G(u), 

where the functional T quantifies in some way the distance between 
the two probability measures \i and v according to a transport cost 
criterium, and F and G are functionals over the space V(Q) (the space 
of probability measures over a domain Q) with opposite behaviour: 
the first prefers spread measures and penalizes concentration while the 
latter, on the other hand, prefers concentrated measures. Obviously 
there are lots of sub-problems which may be of interest, for instance 
the minimizations of the two separate functionals 

$ v (jji):=T(jji,v) + F(ii) and = T(p,,v) + G(v), 

where each time one of the variable is frozen. Also imposing constraints 
like F(n) < H, G(v) < L . . . instead of adding penalizations in the 
functionals may sometimes be considered (and this is in fact the same 
as adding penalizations through some 0/ + oo functions). 

These minimization problems are likely to appear in several phenom- 
ena both in nature and in decision science. For instance in jl] and 
[Hj these variational problems have been intended for urban planning 
models, where the spread measure /x stands for inhabitants, the con- 
centrated v for services and they have to be close in a transportation 
distance sense. Recently a possible choice of the functional §"„ has been 
proposed as a model for the formation of a certain kind of leaves: if 
v = 5q representes the point where feeding arrives to the leaf, then the 
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shape of the leaf is such that it optimizes the quantity of light it re- 
ceives from the sun taking also into account the effective trasport cost 
for bringing feeding all over its shape. 

We present here some choices for the functionals T and G. The choice 
of F is in fact easier since a very good class of concentration-penalizing 
functionals is given by local convex functionals over measures, for in- 
stance 



for any convex fucntion / with /(0) = which is superlinear at infinity 
For these functionals we refer to pp. Here £ is a reference measure 
that may be chosen as the Lebesgue measure C d if Q C M d . By Jensen 
inequality, spread measures with constant density are optimal for these 
functionals. 

Possible choices of T are the following: 

• terms involving Monge-Kantorovich optimal transport cost, as 
in Wasserstein distances T(fj,, v) = W^(fJ,,u) (Wasserstein); 

• terms taking into account traffic congestion effects in trasporta- 
tion distances, as in |4 where the choice T(/x, v) = \\fi—v\Wi for 
a vector space X C H 1 ^) is detailed and justified (congestion); 

• terms reflecting the natural ramified structure of a transporta- 
tion network as in jHj, where a new distance on probability 
measures is introduced according to this criterium (branching) . 

This last possibility is the most suitable for model involving natural 
branching structures like leaves, while the first two seem to be quite 
natural in urban planning. 

For the functional G, before presenting a list of examples, we give a 
possible definition of the concept of concentration-preferring: 

Definition 1. We say that G : V(Q) — > K is a concentration-preferring 
functional if it holds G(t$v) < G(v) for any measure v G V(Q) and any 
1— Lipschitz continuous map t : Q — > Q. 

It is easy to show that any G with this property is minimized by any 
measure 5 X0 , with xq G Q. We list here some funcionals satisfying this 
definition: 

• G(v) = (jspt(V) (atomic), as in location problems, where the 
corresponding T is usually of Wasserstein type; 

• (subadditive, see [T]) 







otherwise 




for a subadditive function g with g(0) = and g'(0) = +oo; 

2 



• G{y) = inf {?-^ 1 (E) |spt(z/) cE,S compact and connected } as 
in the irrigation problem, [2], where T = W\ and a constraint 
on G is considered {length); 

• G(u) = f n Q w(d(x,y))v(dx)v(dy) for an increasing function w 
(interaction) . 

Actually the two first terms are decreasing under the effect of any 
map t and not only under 1— Lipschitz ones. The last choice for G is 
in fact a well-known functional on probability measures which repre- 
sentes interaction energy (or costs) between the particles composing 
v. It has been first studied in a transportation framework by McCann 
in|H], where displacement convexity results are given, with the aim of 
producing uniqueness for variational problems. 

In |3j and [I] two combinations of these functional have been studied 
for urban planning purposes: the Wasserstein + subadditive and the 
congestion + interaction cases, respectively. The congestion + subad- 
ditive case has been exluded in j3] since it leads to a trivially infinite 
functional, and so in this paper we analyze the remaining one, i.e. the 
Wasserstein + interaction case. Many ideas are taken from jl], up to 
the fact that elliptic regularity is replaced by considerations on Monge- 
Ampere equation. Moreover some extra devices are performed and a 
careful use of Monge-Kantorovich theory is needed. 

2. Preliminaries on Optimal Transportation 

We recall in this section the tools and definitions we need in the 
sequel about optimal transport theory and Monge-Ampere equation. 
Our main reference is [Zj. 

Definition 2. Given two probability measures \x and v on a space Q 
and a l.s.c. cost function c : Q x Q — > [0, +oo] we consider the problem 



and the minimizers for this problem are called optimal transport plan 
between \x and v. Should 7 be of the form (id x £)p/x for a measurable 
map t : Q — > Q, the map t would be called optimal transport map from 
fi to v. 

An important tool will be duality theory and to introduce it we need 
in particular the notion of c— transform (a kind of generalization af the 
well-known Legendre transform). 

Definition 3. Given a function x : — > IR we define its c— transform 
(or c— conjugate function) by 




(2.1) 



X c (y) = mi c(x,y) - x{x). 
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Moreover, we say that a function if) is c— concave if there exists \ with 
if) = X c an d we denote by \& c (f2) the set of c— concave functions. 

It is well-known a duality result stating the following equality: 

mm(K) = sup ipd/j + ip c du. (2.2) 
V>e*c(n) Jn Jn 

Definition 4. The functions if) realizing the maximum in (|2.2jl are 
called Kantorovich potentials for the transport from fj, to v. 

Since we will use c(x,y) = \\x — y\ 2 , let us denote by \T/ 2 (^) the set 
of c— concave functions with respect to this quadratic cost. It is not 
difficult to check that 

if) G ^2^) ^ x i— > — — ^ is a convex function on R d restricted to £1. 

Notice that on a bouded Q with diameter D any if) e 1 i r 2(^) is in fact 
2D— Lipschitz continuous. We summarize here some useful results for 
the quadratic case c(x,y) = \\x — y\ 2 . 

Theorem 2.1. Given fi and v probability measures on a connected 
Q C M. d there exists unique an optimal transport plan 7 and it is of the 
form (id x t)$n, provided // is absolutely continuous. Moreover there 
exists also at least a Kantorovich potential if), and the gradient Vif> is 
uniquely determined fi—a.e. (in particular if) is unique up to additive 
constants, provided the density of n is positive a. e. onVL). The optimal 
map t and the potential if) are linked by t(x) = x — Vif)(x) and so t is 
the gradient of a convex function. Moreover it holds if>(x) +if c (t(x)) = 
c(x,t(x)) for fi—a.e. x. 

Starting from the values of the problem (K) in (|2.1|) we can define 
a set of distances over V(VL). For any p > 1 we can define 

W p (fjL, v) = ( min(K) with c(x, y) = \x — y\ p ) 1 ^ P , 

and obviously we will restrict our analysis to the case p = 2. We recall 
that it holds, by duality formula, 

~Wl(n,v)= sup / ipdv+ I ip c dfi. (2.3) 

Theorem 2.2. If Q is compact, for any p > 1 the function W p is 
in fact a distance over V(Q) and the convergence with respect to this 
distance is equivalent to the weak convergence of probability measures. 
In particular any functional \x 1— > W p (/i, v) is continuous with respect 
to weak topology. 

The next step of our analysis is concerned with some regularity prop- 
erties of t and if) (the optimal transport map and the Kantorovich po- 
tential, respectively) and their relations with the densities of \x and v. 
It is easy, just by a change-of- variables formula, to transfrom, in the 
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case of regular maps and densities, the equality v = tj/i into the PDE 
v(t(x)) = u(x)/\Jt\(x), where u and v are the densities of /i and v and 
J denotes the determinant of the Jacobian matrix. Recalling t = V0 
with convex, we get the Monge- Ampere equation 

where M denotes the determinant of the Hessian M0 = det Hcf) = 
det d ®Q x . ■ This equation up to now is satisfied by = id — if> just 
in a formal way We define various notion of solutions for ()2.4)1 : 

• we say that satisfies (|2.4j) in the Brenier sense if (V0)jM- C d = 
v-C d (and this is actually the sense to be given to this equation); 

• we say that satisfies (12 .4|) in the Aleksandrov sense if Hcf), 
which is always a positive measure for convex, is absolutely 
continuous and its density satisfies ()2.4|) a.e.; 

• we say that satisfies ()2.4|) in the viscosity sense if it satisfies 
the usual comparison properties required by viscosity theory 
but restricting the comparisons to regular convex test functions 
(since M is in fact monotone just when restricted to positively 
definite matrices); 

• we say that satisfies ()2.4j) in the classical sense if is C 2 and 
the equation holds pointwise. 

Notice that any notion but the first may be also applied to the equation 
M0 = /, while the first one just applies to this specific transportation 
case. The results we want to use are well summarized in Theorem 50 
in|Zj: 

Theorem 2.3. If u andv are C°' a (Q) and are both bounded above and 
below on the whole Q by positive constants and Q is a convex open set, 
then for the unique Brenier solution of (12. 4j) it holds G C 2 ' a (Q) D 
C 1,a (f2) and satisfies the equation in the classical sense (hence also 
in the Aleksandroff and viscosity senses). 

3. Optimality Conditions for the Interaction Case 

We are here concerned with the minimization problem for the func- 
tional when the transport term is given by T(/x, v) = ^W 2 (fi, v) 
and the concentration one is an interaction term of the form 

G{u)= f V(\x-y\ 2 )u(dx)u(dy), (3.1) 

JSJxSl 

with V : [0, +oo[^ [0, +oo[ a regular increasing fuction. From now on 
Q will be the closure of a convex non-empty open set in M. d . 

A priori a minimizer for this functional may be an arbitrary proba- 
bility measure on the set Q, even a singular one. Our goal is to prove, 
under suitable assumptions and by means of optimality conditions and 
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of an approximation process, that it is in fact an absolutely continuous 
measure with bounded density. 

We provide here an easy optimality condition for the minimization 
of We do not go into details in the computation both beacuse it 
is not difficult to perform it and because it follows the same, trivial, 
scheme as in [H] . Notice that the approximation argument by measures 
with positive densities that we are going to use works in this case too, 
while the alternative proof by convex analysis that may be found in [Hj 
does not, simply because there is no convexity in the term G. 

Theorem 3.1. If a probability measure v £ V(Q) is a minimizer for 
3^, then there exists a constant m such that 

if) + T v > m; ip + T u = m u-a.e., 

where ip is a Kantorovich potential for the transport from v to fx and 
we define 

T u (x)= [ 2V(\x-y\ 2 )v(dy). 
Jn 

Proof. Let us start from the case when fx is absolutely continuous with 
positive density. In this case we perturbate an optimal measure v into 
v t — v + t{i>i — v) for an arbitrary v\ £ V(fl): if we call if>t the unique 
Kantorovich potential from v t to fx which vanishes at a certain fixed 
point xq £ Q, we get (by means of duality formula) 

4>tdv t + / i>tdfx+ / V(\x - y\ 2 )is t (dx)is t (dy) 



> / ip t dv+ / il>%dti+ / V(\x -y\ 2 )u(dx)u(dy). 
Jq Jn Jnxn 

After erasing the term J n ip^ d\x and dividing by t we pass to the limit, 
and we know that ip t converges uniformly (by Ascoli-Arzela) to the 
unique Kantorovich potential ijj from u to fx vanishing at Xq. This 
provides, at the limit, 



/ (V> + T v ) dv x > I {ij; + T v ) dv. 
Jn Jn 



Being V\ arbitrary we get that v— a.e. the function if) + T v must be 
equal to its infimum, and this is the thesis. 

To generalize the result to an arbitrary measure fx, just proceed by 
approximation. This can be performed standardly as in f3] and pro- 
vides the same formula where if> becomes just one of the possibly many 
Kantorovich potentials instead of the only one. □ 

The problem in this condition lies in the fact that the measure v 
appears only in a very implicit way (both in if) and in T v ), and this 
does not allow to derive any estimate on it. We will consequently need 
to pass through an aproximation process, exactly as in jlj. Fixed a 

6 



minimizer v for we will consider a sequence of problems (P e ) e given 
by the minimization of 

V(0) T(fi £ , v) + G{y) + S e A(u) + eW%(v, v e ), 

where 

• (// e )e is a sequence of probability measures approximating // 
with Lipschitz continuous strictly positive densities u £ ; 

• the functional A is given by 

A(u) = l^ a{v)dCd itv = v - CA ' 
[ +00 otherwise, 

for a convex function a : [0, +00 [— > [0, +00] which is superlinear 
at infinity and blowing up at 0, i.e. lim^ + a>(t) = +00, but 
finite and C 2 with a" > c > on ]0, +oo[ (for instance a(t) = 
t 2 + l/t); 

• (S E ) e is a suitably chosen sequence of small positive numbers; 

• {p £ ) £ is a suitably chosen sequence of measures approximating 
v. 

We will prove, exactly as in that this sequence of problems ad- 
mits an uniform L°° bound for their solutions and that we can choose 
the parameter so that these solutions converge to v, thus obtaining an 
L°° estimate for v. The existence of solutions for P £ is trivial by the 
semicontinuity of each term in the sum with respect to the weak con- 
vergence of probability measures on the compact set Q (and moreover 
any term but A is actually continuous, while A is semi continuous by 
convexity) . 

Lemma 3.2. Suppose that \x is absolutely continuous with an L°° den- 
sity bounded by M and that V is a C 2 function with V > 0, then any 
solution v £ to the problem P £ is absolutely continuous and its density 
is bounded by a constant C depending only on M, d and V . 

Proof. First we notice that, thanks to the presence of the term A(u) 
in the minimization problem, v £ must be absolutely continuous with 
strictly positive density almost everywhere. We write down the opti- 
mality conditions for v £ with respect to variations of the form v t = 
v £ + t{y\ — u £ ). From easy computations we get 

ipe + T Ue + 5 £ a'(u £ ) + exs = m £ q.o., 

where ip £ is the Kantorovich potential for the tansport from v £ to fi £ 
and Xe from v £ to v £ (they are unique up to additive constants) and 
m £ is a suitable constant. We get equality almost everywhere due 
to the fact that we already know that v £ > (we identify measures 
and their densities in this context). Since Kantorovich potentials are 
Lipschitz functions and T Ve shares the same regulartiy of the integrand 
(x, y) 1— »• V(\x — y\ 2 ), which is C 2 and then Lipschitz, we get that even 
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a'(u e ) is Lipschitz continuous, and in particular it is bounded. This 
prevents v e to be close to since it holds lim f ^ + a 'if) — —oo. Thus 
we get v e > c £ > 0. Moreover, a!(v s ) is Lipschitz continuous and, being 
a" strictly greater than a positive constant, also the inverse of a' is 
Lipschitz. This proves that u e is a Lipschitz continuous function. We 
can now use regularity theory on Monge- Ampere equation to get ip G 
C 2 ' a (Q) fl C 1,a (Q), since both v e and /i e are bounded both from above 
and from below by positive finite constants (depending on e, anyway) 
and are Lipschitz continuous. The same is true for the Kantorovich 
potential Xe by replacing \i e by v £ . What we can do now is looking for 
a maximum point x of v £ . Notice that such a point will be a minimum 
point for ip £ + T Ue + ex £ - First we prove that x dQ. To prove this 
it is sufficient to prove that the gradient of ip £ + T„ e + ex £ is directed 
outwards at any point of dfl, i.e. V(^ £ + T Ve + ex £ )(xo) • n(xo) > for 
any xq G dQ, where n is the outward normal vector. From the fact that 
the optimal transport map t from v e to fi £ is given by t(x) = x — Vip{x) 
we know that x — Vip(x) G Q for almost any x G Q. In this case, due 
to continuity up to the boundary of V^, this holds for any x and also 
for xq G dfl and implies V^fxo) • n(x ) > 0. Analogously we get 
Vx(x ) • n(x ) > 0. For the gradient of T Ve it holds 

VT„ £ (x )= / W(\x -y\ 2 )(x -y)v £ (dy), 
Jn 

and so VT Ue (xo) ■ n(xo) > since V > and for almost any y G Q it 
holds (a;o — y) ■ n(xo) > 0. This proves that xq lies in the interior of Q 
and this allows us to look at the second derivatives. Taking Hessians 
we have 

H^ £ (x ) + HT Ve {x Q ) + eHxe(x ) > 0, 

where the letter H denotes Hessians and the inequality is in the sense 
of positive definite symmetric matrices. Thus we get 

H^ E (x )>-I (2||V|| c2(n) + £ ), 

since the second derivatives of T Ue may be estimated by those of V 
and from the fact that x 2 /2 — x( x ) is convex we deduce Hx < I- 
This is a uniform estimate from below for Hip £ (xo) and for the con- 
vex function given by <p(x) = x 2 /2 — ip £ (x) we have H(J)(xq) < 
7(1 +£ + 211^11^(0)) • This implies M<P(x ) < (1 + e + 2| \V\ \ c ^n)) d , 
and, from v £ = // £ (V0)M0, we get, for e < 1, 

max v £ = u £ (x ) < 2 d M (l + | \V\ \c^n)Y , 
which is the desired estimate. □ 

Remark 1 . The simple proof we give here of the fact that the gradient 
is directed outwards and no maximum point is allowed on the boundary 
could be used similarly in jlj, thus getting rid of the strict convexity 
assumption in Theorem 6.5 and of the heavy proof of Lemma 6.6. 
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Notice that it could be possible to get the same result even without 
C l regularity for the potentials, just making the proof a bit trickier. It 
would be sufficient to evaluate the increments of the potential in small 
balls around Xq where the gradient is almost everywhere defined and 
such that x — Vip(x), x — V%(a;) G Q for a.e. x. 

Lemma 3.3. It is possible to choose the parameters for the problem P £ , 
i.e. the numbers S £ and the measures v £ and /x e so that any sequence 
of minimizers (u e ) e converges to v. 

Proof. It is sufficient to choose for v £ a sequence of absolutely continu- 
ous measures with Lipschitz continuous strictly positive densities such 
that F(p £ ) < F(u)+e 2 . Then we have A(u £ ) < +00 and we may choose 
5 £ = e 2 A{u £ )~ l . For (/x e ) e we can chose any sequence of absolutely con- 
tinuous measures with Lipschitz continuous strictly positive densities 
approximating /x in such a way that W2(n £ ,fi) < s 2 . Then we have 

T(ji s ,v s )+G{v s )+5 e A(v e )+eW*{y e ,v e ) < T(fi e , T> £ ) +G(v £ ) +5 £ A(v E ) , 

which implies 

F(u £ )+5 £ A(u £ )+eW 2 (u £ ,u £ ) < F{u £ ) + 4L>W 2 (/x £ , /x) + S e A{u e ) 

< F(v)+e 2 + ADe 2 + e 2 

< F(v £ )+e 2 (2 + 4D). 

In the end this implies ^(^o^) < C-Je and, since v £ — ^ P, we get 
v £ ^v. □ 

We can now state our main result and its consequence in the mini- 
mization for the whole functional 5- 

Theorem 3.4. Given a compact convex set Q C M. d with nonempty 
interior and a probability measure /x G L°°(Q), if the function V ap- 
pearing in the definition of the functional G is C 2 and V > 0, then the 
minimization problem for the functional ^ over the space V(Q) ad- 
mits at least a solution and any minimizer belongs in fact to the space 
L°°(Q). 

Proof. As usual the existence is trivial due to continuity and compact- 
ness of V{VL) while, for the L°° regularity, just apply Lemma l3~2l and 
Lemma IH.HI □ 

Corollary 3.5. Given a compact convex set Q C IR d with nonempty 
interior, a C 1 strictly convex and superlinear function f and a C 2 
function V with V > 0, then the minimization problem over the space 
V(VL) 2 for the functional ^ip,u) = |W 2 2 (/x, v) + F{p) + G{y), where 
F is defined by (jl.lj) and G by ()3.1j) admits solutions and, in any 
minimizing pair (/x, u), both /x and v are in fact absolutely continuous 
measures /x = u • C d , v = v ■ C d , with u G C°(Q) and v G L°°(Q). 

9 



Proof. After the usual proof of existence by the direct method in Cal- 
culus of Variations, we refer to [3] for the regularity results on //. Since 
such a measure turns out to be absolutely continuous with continu- 
ous density (hence bounded), we may apply Theorem 13.41 to get the 
regularity on v. □ 

4. An Explicit Example 

In this section we come back to the whole problem of minimizing 
in a very particular case, where we can provide almost explicit densities 
fo the solutions. We consider the case 

. T(ji,u) = \W%M and G{u) = J Qxn V(\x - y\ 2 )u(dx)u(dy), 
as in the previous Section; 

• V(\x — y\ 2 ) = k\x — y\ 2 and so, setting bar(z/) = J n yu(dy), we 
have T u (x) = X\x\ 2 — 2Xx ■ bar(V) + A J Q \y\ 2 i>(dy); 

• F(fi) = |||A*||x,2(n), a particular case of what considered in [Sj. 

The framework we obtain is very similar to the one in 

Theorem 4.1. In the specific case described above, any pair of raini- 
mizers (/i, v) is shaped as follows: 

• fi is concentrated on a ball B{xq, r) (intersected with Q) and has 
a density u given by u(x) = 2^+r( r2 — \ x ~ x o\ 2 )> 

• v is concentrated on the ball B(xq, rj (2A+1)) and it is the image 
of /i under the omothety of center xq and ratio (2A + l) -1 ; 

• Xq is the barycenter of v. 

Proof. First we write down the optimalty conditions given by Theorem 
13. II for the minimization in v with fixed /i and by [3 for the minimiza- 
tion in \i for fixed v. We denote by u and v the densities of \i and u, 
respectively. We may suppose that the barycenter of v is the origin, 
thus obtaining T v (x) = \\x\ 2 + c. We have 

u(x) + ip(x) = C\ a.e. on u > 0; 
if)(x) + \x 2 = C2 a.e. on v > 0. 

Here ip and if) are Kantorovich potentials for the transport from [i to v 
and from v to //, respectively. From the second condition we can infer 
Vip(x) = —2Xx a.e. on v > 0. Being v absolutely continuous, this 
equality is valid v— a.e.. This means that the optimal transport map t 
from v to n is given by t(x) = x — Vip(x) = (2A + l)x. By uniqueness 
of the optimal transport plan, which is in this case expressed both as 
a transport map from v to fi and viceversa, we know also the optimal 
transport map s from \i to v which will be s(x) = x/(2X + 1). From 
duality theory in optimal transportation we know the following equality 

<f(x) + ip c (s(x)) = c(x, s(x)) = -\x - s(x)\ 2 , 
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and thus we get u(x) — c\ — \\x — s(x)\ 2 + ip c (s(x)). Since ip c is a 
Kantorovich potential from v to /i, we know that it must agree (up 
to constants) with ip on any connected component of {v > 0} (notice 
that in this case we are allowed to speak about connected components 
of this set, since {u > 0} is an open set, being u continuous, and the 
set {v > 0} is just an omothety of it). So, let uj C VL be a connected 
component of {u > 0}. On (2A + l) _1 u; we have tp c = ip + c 3 and so we 
get 

u{x) = C4 — -\x — s(x)\ 2 + ip(s(x)) = c 5 — \x\ 2 



u(x) 



c — IX 



2 



(4.1) 



2' v " rv v " 0,1 2A + 1 
From this expression it is clear that duj \ dQ (where u must vanish) is 
contained in a circle line around 0. This implies that belongs in fact 
to u, since no boundary of u is allowed in the interior of a certain ball 
around 0. So there is in fact just one connected component for {u > 0} 
and so we get 

A ' 
2A"TT. 

From this it is easy to recover the density of v since v = sj/i and we 
get the thesis. The point xq wich turns out to be the center of the balls 
which are supports for fi and v is in this notation 0, the barycenter of 
z/, as in the thesis. It is clear that in this case \i and v share the same 
barycenter since they are omothetical. □ 

Remark 2. In this example the density v shares the same regularity of 
u except at the points corresponding to boundary points of Q where u 
is positive, i.e. if at x G dQ it happens u(x ) > then at s(x ) we have 
a jump for v. It is clear from the fact that u is 2A/ (2A + 1)— Lipschitz 
continuous (it follows from the explicit formula) that we have, recalling 
also j n u dC d = 1, 

1£ ( i " fM+ 2fTT c ) 

where D is the diameter of Q. This implies, for small Q, inf u > 0. In 
this case u would be positive at any point of dfl and v discontinuous 
at any point of s(dQ). This gives examples when the L°° regularity for 
v is quite sharp (in the sense v C°(Q)). 

Remark 3. In the explicit example above there remain to be deter- 
mined both the constant c in the formula for u and the position of the 
barycenter xq. In some simple cases this too is possible. Notice that, 
once fixed Xo, the constant c may always be recovered by imposing the 
condition of being probability measures. For instance if Q is a ball, 
we may see that xq may not be the barycenter of a density u shaped 
as in (j4.1|) unless the set B(xo, 2A _1 (2A + 1)) D ft is a ball around xq. 
This happens for large Q whenever the ball _B(xo, 2A _1 (2A + 1)) does 
not touch the boundary dQ or, in general, when xq is the center of the 
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ball Q. In the first case (Q a large ball) we have several solutions for 
the problem (non-uniqueness), obtained from each other under trans- 
lations, and u and v are continuous; in the second (Q a small ball) we 
have uniqueness of the solution, with u a radial continuous function 
around the center and v a rescaling of u on a smaller ball. 
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